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Abstract. Given a graph Laplacian with positively and negatively weighted edges we are inter¬ 
ested in characterizing the set of weights that give a particular spectral index, i.e. give a prescribed 
number of positive, zero, and negative eigenvalues. One of the main results of this paper is that 
the set of signed Laplacians that exhibit multiple zero eigenvalues is “small”, and that eigenvalue 
crossings are nongeneric — specifically, eigenvalues repel each other near zero in a sense that can be 
made precise. We exhibit an algebraic discriminant that measures the level of repulsion, and show 
that this discriminant admits a combinatorial interpretation. Conversely, we exhibit a constructive 
method for finding the sets of Laplacians that exhibit a large degree of degeneracy (many eigenvalues 
at or near zero) in terms of these discriminants. 
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1. Introduction. Let T = { 7 ^} be a symmetric weighted graph, let £(r) be 
its graph Laplacian, and let n+(r) be the number of positive eigenvalues of this 
Laplacian. It is a classical result that if 7 ^ > 0, then ro+(T) = 0. In I], the authors 
computed bounds for n+(T) when 7 ^ are allowed to be both positive and negative. 
We showed that the sign topology of the graph (i.e., only knowing which edges were 
positive and which were negative) determines strict upper and lower bounds on n + (r), 
and, moreover, that these bounds are saturated: there is some choice of weights that 
achieves these bounds. 

There are graphs for which n+(r) is independent of the magnitude of all of the 
weights (these graphs were termed “rigid” in |l]); however, for “ most” graphs, n+(r) 
depends on the magnitudes of the weights attached to each edge. A natural question 
is then, for any particular sign topology, to describe the set of weights which give a 
prescribed number of positive eigenvalues. In this paper, we study this question in 
its most general formulation, describing these sets for any particular sign topology, as 
well as considering questions of genericity. 

We also consider applications of this question to designing networks with partic¬ 
ular dynamical properties, in particular those that can differentiate multiple signals, 
and those that are stable but can support multiple solutions on long timescales. 

1.1. Results of paper. If V = { 7 ^} be a weighted symmetric graph, then we 
define T(L) to be its graph Laplacian, i.e. £(r) is the matrix whose entries are 


£(i 'h 


’liji *7 ~ ji 

~ L/c i=j- 


( 1 . 1 ) 


We define (n_(r),no(r),n_|_(r)) to be the number of negative, zero, and positive 
eigenvalues of £(r). 

Let G = (V,E) be a symmetric unweighted graph, where all of the edges are 
colored black or red. We term this colored structure the “topology” of the graph 
G. We also define the two subgraphs G + (resp. G_) to be the subgraphs where we 
consider only black (resp. red) edges. We define c(-) of any graph to be the number 
of its connected components. 
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Let us now choose the convention that black edges will correspond to positive 
weights, and red edges to negative weight^ It is clear that, given a weighted graph T, 
there is a unique colored graph G associated with it — we simply forget the magnitude 
of the weights and retain only their signs. Conversely, for any signed graph G, there 
is a natural map from the positive orthant to the set of weighted graphs T. 

We then show a number of results about the index of all weighted graphs with 
a given topology. We first state the following results, which are a generalization and 
slight restatement of the authors’ previous results in 11: 

1. If we fix the positive weights of the graph and choose the negative weights 
in a neighborhood of zero, then generically the Laplacian graph has c(G + ) — 
1 positive eigenvalues, a single zero eigenvalue, and the rest positive. In 
particular, if G+ is connected, then for negative weights sufficiently small, 
the Laplacian is negative semidefinite. 

2. If we fix the positive weights of the graph and choose the negative weights in 
a neighborhood of oo, then the Laplacian graph has n —c(G_) positive eigen¬ 
values, a single zero eigenvalue, and the remainder negative. In particular, if 
there are any red edges, then we can destabilize the matrix by choosing their 
weights large enough. 

If we consider a one-parameter family of Laplacians with fixed black weights and 
move the red weights from zero to infinity along some ray, then it follows from the 
above that r := n — c(G + ) — c(G_) + 1 eigenvalues move through zero from left to 
right. The next question is: does this one-parameter family have r distinct crossings 
of individual eigenvalues, or do multiple eigenvalues cross at the same time? (Note 
here that the Laplacian will always have an eigenvalue fixed at 0, so a simple crossing 
corresponds to an eigenvalue with multiplicity exactly two.) 

To state the answer precisely, we need some notation. Let G be a fixed graph with 
B black edges and R red edges, and denote the weights on these edges by £ R B 
and w _ £ M. R . This induces a map from x to all weighted graphs with a given 
topology, by setting the positive weights to be the values of w + and the negative 
weights to be —W-. If we restrict w+,w- to the positive orthant, then there will 
be positive weights on the black edges and negative weights on the red edges. The 
mapping from x to weighted graphs induces a topology on weighted graphs, 
and it is in this sense that we use the term generic below. 

The one-parameter family of graphs described above corresponds to fixed w+ and 
the ray tw _, t £ [0,oo) — this gives a one-parameter family of graphs, which gives a 
one-parameter family of sets of eigenvalues, indexed by t. 

The main results of this paper are as follows. We will only state and prove the 
results in the case where T is connected (if not, the Laplacian is the direct sum of the 
connected copies, so these results generalize in a straightforward manner). 

1. Fix any u> + for the positive weights. For a generic set of W-, all of the 
eigenvalues have distinct crossings, i.e. there are r distinct values of t at 
which the Laplacian has a single eigenvalue crossing zero. 

2. There exists a generic set of w+ such that for any u>_, all of the eigenvalues 
have distinct crossings. Moreover, for any fixed uj+, if we consider the set 
of all W- with ||u!_|| = 1, then there is a “minimum distance” in t between 
successive eigenvalue crossings. This shows the phenomenon of level repulsion 
for generic Laplacians. 


1 We are stealing this convention from the accounting industry. 
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3. To obtain multiple eigenvalues crossing simultaneously, then from the previ¬ 
ous two statements, one must be “doubly nongeneric” in a certain sense. We 
introduce an algebraic expression, related to the discriminant of a polynomial, 
that determines when this nongeneric situation can occur. This discriminant 
can be expressed in terms of the homology of the graph, and corresponds to 
a signed count of certain spanning 2-forests in the graph. 

4. Finally, we give sufficient f? 1 conditions on the vector that guarantee 
stability of the Laplacian, and moreover characterize these conditions in a 
combinatorial manner. 

There is an obvious duality in the statements above, since the spectrum of — L is 
just the spectrum of L times — 1, so we could also choose to fix the negative weights 
W- and vary w+ and obtain obvious analogues of the above results. 

1 .2. Applications. We present two applications motivating these questions: 

1] Consider the problem of designing a linear network to differentiate multiple signals, 
where the network structure is prescribed. More specifically, consider the system 

(] T 

— =Ax + f(t), x € r, f:R^R n (1.2) 

at 

where /(f) represents an external signal, the matrix A represents the network connec¬ 
tivity, and the vector x(t) represents the network response. In order for the system 
to be stable we require that the spectrum of the matrix A lie in the left half-plane. 
It is clear that eigenvectors corresponding to eigenvalues far from zero are difficult to 
excite, so the main response from such a system is from the eigenvalues near zero. 
For such modes the system above acts as an integrator, integrating the projection of 
the external signal onto the eigenmodes. Being able to recognize multiple signals is 
the same as saying that we would like to choose A to have a large dimensional kernel 
one dimension for each signal we would like to be able to recognize — or, perhaps, 
to choose A so that it has many eigenvalues near zero so that the responses track the 
signal with a slow decay. 

Of course, it is not difficult to design a linear system with spectrum wherever 
we would like: simply choose the eigenvalues, then any matrix similar to the corre¬ 
sponding diagonal matrix would do. However, notice that in general this gives a dense 
matrix, and naively it is not clear how one can choose the eigenstructure so that the 
eventual linear system is compatible with a desired topological structure. However, 
using the results of this paper, we show how this can be accomplished: if we can 
design a network with a fc-fold degeneracy at zero, then an open set of perturbations 
of the weights of such a system will give the desired network. 

2] Given a graph G = (V. E) and symmetric coupling functions = tPjii'), define 


—Xi = Fi(x) := Ui + ^ Vijixj-Xi). 


(1.3) 


One famous case of this model is the Kuramoto oscillator network [2jj5], where = 
7 ij sin(-). Assume x* is a fixed point for (1.3), i.e. Ti(x*) = 0 for all i. The stability 


of this point is determined the index of the Jacobian J, where 


<Pij( x j -x *)> 

“E fc¥4K - 
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In the Kuramoto case, the off-diagonal terms are given by 7 * j cos (at* — a:*). The 
Jacobian J is a graph Laplacian of the form (1.1); thus, determining the stability 
indices for fixed points of (1.3) is related to the problem studied here [6-10 . Of 
course, identifying those fixed points whose Jacobian is negative semi-definite gives 
the attracting fixed points for the system, but in fact being able to determine which 
of these points have one unstable eigenvalue is important to understand metastable 
transitions for stochastic versions of this system 11 . 

If the components of x* are close enough to each other, then cos(x* — x*) > 0, 
and the Jacobian is negative semidefinite by the classical theory (see Theorem 3.1 
of [l2]). However, one might ask about the stability of “splay states”, i.e. those 
stationary points where some of the components are far enough from each other to 
make the cos(-) term negative. For a generic choice of uj 1 , ip[j{x* — x*) is non-zero 
for all (i,j) £ E, implying that the graph determining J and the graph defined by 
the original interactions in (1.3) have the same underlying topology. Thus we have a 
fixed network topology, and want to understand the effect of some edge weights being 
negative. The boundary of the region where the matrix is negative-semidefinite with 
a one-dimensional kernel is, of course, the set of points where the matrix is negative 
semi-definite with a higher dimensional kernel 
motivation for studying this problem. 


10 , and this provides yet another 


2. Statement of main results. In this section, we present the main results of 
this paper, leaving the proofs for later sections. Many of the definitions in this section 
are identical to those of 1 , but we include them here for completeness. 

2.1. Weighted graphs, signed graphs, and the Laplacian. Definition 

2 . 1 . 

• A graph G = (F, E) is a set V of vertices and a set E C F x F of edges. 

• A signed graph is the triple G = (F, E, a) where (F, E) is a standard graph, 
and with a map <r: E —> { red , black}. 

• A weighted graph is the pair T = (F, { 7 ^}) where 7 ^ £ K. The edges o/T 
are those ( i,j ) with 7 \j 7 ^ 0 , and we say that 7 ^ is the weight of edge i ++ j. 

Any weighted graph corresponds to a signed graph in an obvious manner; we will 
colloquially call this the “topology” of the weighted graph. 

Definition 2.2. Given a weighted graph T, the Laplacian o/T is the matrix 
£(T) with 


C ■ = 


Tiji * 7^ j j 

i = j- 


The index of (the Laplacian of) T is the triple of integers 


(n_(r),n 0 (r),n+(r)) 


( 2 . 1 ) 


( 2 . 2 ) 


giving the number of negative, zero, and positive eigenvalues of £(T). 

Notation 2.3. If G is a signed graph, we will denote by G+ the subgraph con¬ 
taining only the plus edges, and G- the subgraph containing only the minus edges, 
with a similar convention for T± when considering weighted graphs. We also denote 
c(G) as the number of connected components of a graph, so that c(G+) is the number 
of connected components of a signed graph when we only consider positive edges, etc. 
We only consider connected graphs in this paper, so that c(G) = l, but we allow for 
c(G+) and c(G-) to be larger than one. IfT is a weighted graph, then we can associate 
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it to a signed graph in the obvious way, and thus all of the notions above make sense 
as well, i.e F +; c(T + ), etc. 

2 . 2 . Crossing polynomial and spectral variety. Definition 2.4. IfT is a 

tree, define 7 r(T) to be the product over the edge weights in the tree 

AT) := I] TH- (2-3) 

Let r be a weighted graph with |V(T)| = N, and £(T) be its graph Laplacian. We 
know that £(T) has a zero eigenvalue, and therefore det(£(T)) = 0. Order the n 
eigenvalues of C(T) so that Ai = 0, then define 

M(r)= [ ^ ( 2 -4) 

*=2 

Thus A4(F) is proportional to the linear term in the characteristic polynomial of the 
Laplacian. Also, At(T) ^ 0 iff 0 is a simple eigenvalue of C(T). 

With this notation the Kirchhoff matrix tree theorem can be stated as follows: 
Lemma 2.5 (Weighted Matrix Tree Theorem). Let T be a connected, weighted 
graph, and ST(T) the set of all spanning trees ofT. Then 

M(T)= J2 (2-5) 

Tesr(r) 


Remark 2.6. This is Theorem VI.29 in the text of Tutte '131: a proof is provided 
there. Notice that if all of the edge weights are non-negative, then the sum in (2.5) is a 
sum of positive terms. This is an alternate proof that the kernel of a graph Laplacian 
with positive weights is negative semidefmite, and has a simple kernel, when the graph 
is connected. However, once we allow negative weights, the sum on the right-hand 
side can have cancellations and will not be sign-definite. 

We will denote the vector w = ( w + ,w _) £ (M) B x (M)^ by 


(si, S 2 , ■ • ■ , S B , t\, t 2 , ■ - • ,Fr). 


The s k are the weights of the black edges, and the weights on the red edges are — tk- 
Every spanning tree will have N — 1 edges. Choose a spanning tree T and for 
some k, it will have k red edges and N — k — 1 black edges. The contribution to the 
crossing polynomial is a term of the form 

7t(T) = (-l) fe s il s i2 . .. s iN _ k _ 1 t 1 t 2 ■ ■ ■ t k . (2.6) 

If we denote STk as the set of spanning trees with exactly k red edges, then 

N — l 

M(T)= J2 AT)=J2 E (-1 ) k s il s l2 ...s iN _ k _ 1 t 1 t 2 ...t k . (2.7) 

Tesr(r) fc=0 reSTfc(r) 

We see from this that At(T) is a homogeneous polynomial of degree N — l with 
alternating signs. 

For the remainder of this paper, we will think of the sV s as parameters and the 
t k s as variables. Then At(T) is a polynomial in the variables (t±,... ,tjf), where the 
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coefficients are set parametrically by combinations of the Sfc’s. This corresponds to 
fixing the black edges with specific weights, and varying the red edges. 

Proposition 2.7. The polynomial A4(F) has terms of degree k iff 

c(G+) - 1 < k < N-c(G-). (2.8) 


Proof. The proof is the same as the proof of Lemma 2.14 from I], which we 
repeat here. Every spanning tree must have at least c(G + ) — 1 red edges to connect 
the components of G+ together. Conversely, it must have at least c(G_) — 1 black 
edges for the same reason. Since every tree has iV — 1 edges, it has at least c(G+) — 1 
and at most N — c{G-) red edges. This gives the upper and lower bounds. 

To see that the polynomial also contains all of the intermediate terms, we can 
argue in two ways. One argument is to use contraction-deletion repeatedly (see Re¬ 
mark 2.17 of |1|) to see that the coefficients satisfy log-convexity; this implies that all 
of the intermediate terms appear in the polynomial as well. A more concrete argu¬ 
ment: we can note that for any number k with c(G+) — 1 < k < N — c(G_), there 
is a subforest of T_ with exactly k edges. Since T + is connected, this means that we 
can extend this subforest to a spanning tree of T using only edges from T + , and thus 
there is a spanning tree of T with exactly k edges, giving the term of degree k in the 
polynomial. □ 

Proposition 2.8. All eigenvalues are nondecreasing as a function of anyth- 

Proof The most straightforward proof is to use the results of 1 . Consider the 
ray in M. R given by 


if l, t 2 , * • •; t(cx i, ■ • •, olr), 

with oti > 0. This reduces to the framework of 1|, where the crossing polynomial 
was a function of only one independent variable, t. Theorem 2.10 of 1 implies that 
along this ray, n_t-(F) is nondecreasing, and Lemma 2.18 of [I] implies that it strictly 
increases by one every time we cross the set A4(T) = 0. Inside this set, there is always 
a double eigenvalue at zero. 

A self-contained argument is as follows: First note that for any x, 
x t C(T)x=- ^2 7 ij^i-Xj) 2 . 

Clearly, if we increase any tk , then this will decrease one of the 7 y , thus increasing this 
quadratic form. From the Courant minimax theorem, this implies that the eigenvalues 
of £(r) are nondescreasing as a function of tk, since all eigenvalues can be written as 
the maximum of this quadratic form restricted to some subspace. □ 

2.3. Coefficients of polynomial. In Definition |2.4| we noted that £(F) has 
a multiple zero eigenvalue iff A4(r) = 0, and from Proposition |2.8| this means that 
eigenvalue crossings occur on the zero variety Z(T) = |A4(r) = 0} in K^, i.e. the 
connected components of Z(r) c in are exactly the sets where n+(r) is constant. 

Definition 2.9. For any spanning tree T ofT, we define 7 f(T) as the product of 
the weights of all of the black edges in T, i.e. 

e{T)= "Ar 

i<j, 

c i,j)eE(T ), 

7ij>0 
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Definition 2.10. Let r be a weighted graph with R negative edges, and let I,J 
be disjoint subsets of [1?]. We define Sj as the set of all spanning trees of F which 
contain all of the red edges in the index set I, and none of the red edges in index set 
J, i.e. 

Sj (r) := {T | T is spanning tree of T,ei £ T Mi £ /, ej ^ T Vj € J} , 
and we define 

S/(r) := n *(- T). 

Tesf 

Finally, as shorthand, we will define 

Si = Sf, Sj=Sf, 
where the complement I c = [i?] \ /. 

Definition 2.11 (Deletion and contraction). LetT = (V, E) be a weighted graph, 
and e £ E(T) an edge. 

• We denote by F \ e the graph obtained by removing edge e, and call this the 
deletion of edge e from T. 

• If e = (v\,V 2 ) is an edge with v\ ^ V 2 , we define T.e as follows: identify the 
two vertices V\ and V 2 as a single vertex v*; for any vertex w connected to v\ 
or V 2 , we define the new edge weight "f v * tW = Jvi,w + lv 2 ,w- We call this the 

contraction of edge e in T. 

Definition 2.12. Let T be a graph with R red edges. If I, J are two disjoint 
subsets of {1,... ,R}, then the graph Tj is the graph obtained by first deleting all of 
the edges in I, and then contracting all of the edges in J. As above, we write T/ as 
shorthand for Tj . 

Proposition 2.13. Using this notation, 

S I (T) = M(T I ). (2.9) 

The proof is the standard contraction-deletion theorem, q.v. 14j §13.2]. 

Remark 2.14. Let T be a graph with R red edges. For any I C [i?] ; the coefficient 
of the term t 1 , i.e. nif=i <5/(r). In particular, this allows us to use the compact 
notation 


M(r)=Y, s ^ tI - 

ie 2 R 


( 2 . 10 ) 


If we constrain all positive weights to be one, the quantity S/(F) has a the com¬ 
binatorial interpretation as the number of spanning trees that contain all red edges in 
I and do not contain any red edges not in I. 

Remark 2.15. We do not use (2.101 explicitly in the results below, but is quite 
powerful computationally: given T,I, to compute S/(r), we can compute T/, write 
down its graph Laplacian C(Ti), then compute the linear term of its characteristic 
polynomial. This process can be defined solely in terms of matrices, is easily automat¬ 
able using a computer algebra system, and is relatively inexpensive. 

Finally, we state the main theorem of the paper: 

Theorem 2.16. For generic positive weights, all red rays give single eigenvalue 
crossings. For all positive weights, generic red rays give single eigenvalue crossings. 

We prove these theorems in Sections [3] and [4] below. 
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3. Multiple eigenvalue crossings, if R = 2. In this section, we concentrate 
on the case where there are exactly two red edges (see Section [4] for the case of more 
than two red edges). In this section, we write the vector (aq,^) as (x,y). 

Let us assume that r + is connected. Then the main theorem of [lj, Theorem 2.10, 
implies that as the weights on the red edges go to zero, n+(T) = 0, and if they are 
sufficiently large, n+(T) = 2. Only two eigenvalues cross through zero, so either they 
cross as a degenerate pair, or they do not. In this section we compute the conditions 
that determine this. 

If T_|_ is not connected, then the crossing polynomial will not have a constant 


term, but after factoring out the lowest-order terms, it will be in the form (3.11, so 
that analogous statements about genericity also hold. 

3.1. Discriminant. Here we connect the discriminant of the 2-variable crossing 
polynomial to the combinatorics of the graph. Recall from before that we have 


-M(r) = Anxy - A w x - A 01 y + H 00 , 


(3.1) 


where Aj = Sj{T). Since T + is connected, H 0 o > 0. 

Definition 3.1. The discriminant of the polynomial (3.1) is the quantity 


A — AhAqq — Am A 


oiciio- 


Lemma 3.2. If A/0, the zero-set Z{T) defines a hyperbola, with the minimum 
Euclidean distance d between the branches given by 


d = 


V2A 

~A^' 


(3.2) 


In the case A = 0 the set Z(T) is a reducible variety and degenerates to the union 
of the lines x = Am/An and y = A^/An- 

Remark 3.3. This clearly illustrates the level-repulsion phenomenon: while we 
can make one eigenvalue vanish by varying the strength of one of the edges (say x) 
we typically cannot make two eigenvalues vanish by varying the strength of two edges. 
The typical situation, when A^O, is that there are two disconnected branches of the 
zero-set, each corresponding to a single zero eigenvalue. It is only in the degenerate 
case where A = 0 that we are able to create a double zero eigenvalue. 

3.2. Proof of Theorems 12.161 for R = 2. It is clear from these formulae 
that the presence of a potential double root is nongeneric, since the weights of the 
matrix need to satisfy a polynomial equation. Thus for R = 2, we have proved 
Theorem |2.16[ as long as A / 0, then all rays give single eigenvalue crossings, and this 
is clearly generic. Moreover, even if A = 0, the only way to actually obtain the double 
eigenvalue crossing is along the ray that passes through the point {Am/An, A w /An). 

From this, it is clear that a double eigenvalue is doubly nongeneric: we need the 
discriminant to be zero, and even then, we need to choose the correct ratio of weights 
on the red edges to hit the sweet spot. See Figure [3Jj 

3.3. Counting forests. We have proved the main theorems by showing that 
the discriminant is the quantity that controls how closely the eigenvalues can cross 
We now give a combinatorial interpretation to this discriminant. 

Definition 3.4. For any graph, a spanning 2-forest is a spanning subgraph 
with two connected components — in other words, a spanning tree with one edge re- 










n + =1 

CM 

II 

c + 

=0 

+ 3 

II 


Fig. 3.1. A schematic diagram of the nondegenerate and degenerate cases when there are two 
red edges. 


moved. More generally, a spanning fc-forest is a spanning subgraph with k connected 
components. 

Definition 3.5. Let G = (V,E) be a signed graph with \V\ = n. Let U,W be 
subsets of[n ] with \U\ = \ W\ = k. We define as the set of all spanning k-forests 

such that every tree in the forest has exactly one vertex in U and one vertex in W. 

We can associate a sign to each forest in F € F^ w , denoted e(F), as follows. 
Choose and fix some enumeration for U, W, and define the map g: W -A U by saying 
that u = g(w) is the unique element of U that is in the same component of F as w. 
Then g(W) is a permutation of U and thus has a sign, and this is e{F). 

Theorem 3.6. Let T be a weighted graph with two red edges and where the red 
edges have weight x and y. Let i,j be the vertices of the first red edge, and k,l be 
the vertices of the second. We can form the polynomial (3.11 as above. If we choose 
U = {i,j} and W = {k,l}, then 


A = - E < F M F ) 


(3.3) 


. FGF, 


( 2 ) 


where tt(F) is defined as in (2.31. If all of the black edges have integral weight, then 
this and 


( |3.2| ) implies that the Euclidean distance between the two branches of Z(T) 
is y/2 times a rational number. 

To prove this theorem, we first need some auxiliary results. 

Theorem 3.7 (Chaiken 1 15]). Let G = {V, E) be a signed graph with |V| = n. Let 
U,W be subsets of[n] with \U\ = \ W\ = k. Define to be the minor determinant 

formed by removing the rows U and the columns W from the graph Laplacian £(G). 
Then 


\£u,w\ = ^ a(F). 

F ^u% 


(3.4) 
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Proof. See Chaiken 15] §3] for a proof in the general (directed graph) case. He 
uses this result to prove the directed MTT, the directed generalization of (2.5) above. 
□ 


Remark 3.8. The case we are particularly interested in here is \U\ = \W\ = 2. 
In this case, forests which assign the least element of U and the least element of W 
to the same component are counted + 1 , and forests which assign the least element of 
U and the least element of W to different components are counted — 1 . 

Remark 3.9. It is not hard to see that this theorem implies that there are a 
number of linear relationships among the minor determinants of a graph Laplacian. 
One example is 


|-£12,13 1 + |£12,14 1 — |£12,34 1 • ( 3 - 5 ) 

To see this note that the first term, |£i2,i3 1 ? is equal to minus the number of 2-forests 

with vertex 1 in one component and vertices 2 and 3 in the other component. This 

collection of 2-forests can be split into two types: those in which vertex 4 is in the 
component with vertex 1 and those in which vertex 4 is in the component with vertices 
2 and 3. Similarly the second term, |£i 2 ,i 4 1 5 is equal to the number of 2-forests 

with vertex 1 in one component and vertices 2 and 4 in the other component. This 

collection of 2-forests can similarly be split into two types: those in which vertex 3 is 
in the component with vertex 1 and those in which vertex 4 is in the component with 
vertices 2 and 3. The terms in which vertices 2, 3 and 4 share a component cancel, 
and we are left with the number of 2-forests in which 1 and 3 share a component and 
2 and 4 share a component, minus the number in which 1 and 4 share a component 
and 2 and 3 share a component. This is exactly |£i2,34 1 - A similar identity, which we 
will use in the proof, is 

| £l 3,13 1 + | £ 14,14 | + | £l 3,23 | + |£ 14,23 | = | £12,34 1 • (3-6) 

The next result we use is an identity variously attributed to Dodgson 16 , Ja¬ 
cobi 17|, and Desnanot [18| on determinants of minor determinants. There is a nice 
bijective proof due to Zeilberger [19] , 

Theorem 3.10 (Dodgson). Let M be a matrix, and Mu,w denote the submatrix 
formed by deleting the rows in U and the columns in W. Then 


\M\ \M ijM \ - |M i>fc | \Mj t i \ = - \M itl M j>k \ 


This identity underlies Dodgson’s Method of Condensation, an algorithm for ef¬ 
ficient hand computations of determinants. We are now in a position to prove the 
main theorem: 

Proof of Theorem 


3.6, For this proof, we will use the symbol Cjj,w to be the 


determinant of the Laplacian matrix with rows U and columns W removed. 

Let us first assume that the red edges x and y do not share a vertex. Renumber 
the vertices so that edge x connects vertices 1 and 2 and edge y connects vertices 3 
and 4. First note that any 1-minor determinant is (up to sign) P(x,y). In particular 
we have that 


P(x, y) = £ 1,3 = Hoo + A 10 x + A m y + Anxy. (3.7) 


Note that we also have 


£12,13 + £12,23 — A10 + Any 
10 


(3.8) 




and that 


£13,34 + £14,34 — An + A n x 

and finally that 

£123,134 + £124,134 + £123,234 + £124,234 = Ai- 
Some algebra gives that 

Ai(Ao + Ao^T Ai 2/ +Ai^y) — (Ao + Aiy)(Ai + Ai®) = Ai Ao — Ao Ai = A 

or 

A = (£l23,134 + £l24,134 + £l23,234 + £l24,234)£l,3 — (£l2,13 + £l2 ,23) (£l3,34 + £l4,34) 

and applying the Dodgson identity (here £ 1,3 is the matrix M) gives 

A = —£12,34 (£13,13 + £14,13 + £13,23 + £14,23) = — (£l2,34) 2 

In the case where x and y share a vertex we number the shared vertex 1 and the 
other two 2 and 3. Then a similar calculation to the one above gives 

A = — (£12,13) 2 


Corollary 3.1. Assume T has two disjoint red edges, x and y, and denote 
x = ( i,j ) and y = ( k,l ). If there is an automorphism of T that exchanges i,j and 
fixes k,l (or vice versa), then A = 0. 

Proof. This automorphism fixes n(F) but negates e(F), thus giving A = —A. □ 


3.4. Cycles and homology. The expressions derived above give explicit for¬ 
mulae in terms of sums over spanning trees and spanning 2-forests; we now show that 
there is a dual formula in terms of minor determinants of the cycle intersection form 
for the homology of the graph T. 

Theorem 3.11. Suppose that the original graph T has co-rank c = E — n + 1, 
and that removing the edges x and y does not disconnect T. We construct a cycle 
basis for the graph in the following way: let V\,V 2 , ■ ■ ■ w c -2 denote a cycle basis for the 
graph with edges x and y removed, v c -\ denote any cycle through x but not through y 
and v c denote any cycle through y but not through x. 

In other words the basis should be chosen so that there is a unique cycle through 
each of the edges x and y, but the basis is otherwise arbitrary. Next we define the 
matrix F to be the E x c matrix with rows given by {uA=i- £e£ AI) denote the 
minor determinant formed by removing the i th row and the j th column of FF*. Then 

a = im c c - 1 i 2 . 

Proof. The proof here is analogous to a similar result due to Sjogren 20 Theorem 
1] that shows that the determinant of FF l is the number of spanning trees of the 
graph. The proof here is similar in flavor, except that, as we are dealing with a minor 
determinant, the basis must be chosen in the particular way described above. 


Having chosen the basis as above, the main idea is the Cauchy-Binet formula 21 


§0.8.7]. We will have to consider various submatrices of F. Superscripts will denote 
deleted rows of the matrix (at most one row will be deleted), and subscripts will 
denote the retained columns. For instance if 5 is a subset of the edges then FfT 1 will 
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denote the submatrix obtained from F by deleting the (c — l) st row (corresponding 
to the cycle through edge x) and deleting all columns not in S. 


where the sum is over all subsets S of 
Using Lemma 


Then the Cauchy-Binet formula 21 §0.8.7] gives that \Mf x | = 1^5 ||-F 


Sh 


3.12 


the edge set E of size c — 1. 
the remainder of the proof is straightforward: 


IK 


i c— 1 1 


f c q 


unless T/ (S\Jy) and T /(SUx) are both trees. This is clearly equivalent to T / ( SUyUx ) 
being a spanning 2-forest where x and y each have one vertex in each component, so 
| M/r 1 ] reduces to a signed sum over spanning 2-forests with this property, with the 
sign given as in the theorem. □ 

Lemma 3.12. Using the notational conventions established in the statement and 
proof of Theorem \3 .1 1\ we have 

• |FlU | = 0, ±1. The determinant is zero unless T/(SUy) is a tree, in which 
case l-FcU 1 ! = 0, ±1, and similarly for |Fg| = 0 . 

• In the case where E/(S U y) and E/(S U x) are both trees then the sign of 

i n 


|Fg| is determined as follows: E/S is a graph with a unique (up to 


sign) cycle that contains edges 

Y 


and 

TIC—1 


“y” ■ 


If this cycle traverses edges 


F c s ~ L m\ = l, otherwise |F^ i ||F|| = -1. 


Note that this sign is independent of the orientation assigned to edges “x” 
and “y” and the direction of the cycle through them. 

Proof. Let c=E — N+lbe the co-rank of the graph. We first note the following: 
if 5 is a set of c edges of T then \F§\ = ±1 if T/S 1 forms a tree and \F§\ = 0 otherwise. 
This is proven by Sjogren [20j, but for the convenience of the reader we give a quick 
proof here. First, to see that |Fg| = 0 if T/S is not a tree we give an explicit element 
of the kernel. Since T/S is not a tree it must contain a cycle. Express this cycle in 
terms of the basis as Q-iVi. Thus aFg gives the edges of this cycle lying in S. Since 
this cycle lies entirely in the complement, this is zero. 

To see that \Fg\ = ±1 when T/S is a tree T we first consider the case where the 
basis is chosen in a special way: namely for each edge i in S we consider the unique 
cycle in T U i. This gives a cycle basis. In this particular basis the matrix F§ is a 
permutation matrix, since each edge in S lies in a unique cycle. Now any integral basis 
is related to this special basis by a matrix U with \U\ = ±1, and the corresponding 
matrices are related by F§ = UFg. Therefore for any basis \Fg\ = ±1. 

We are interested in the (c — 1) x (c — 1) matrices Fg and FS ~, where S is now 
a subset of c — 1 edges, but these matrices can be realized as minors of a matrix Fg 
of the above form. It is clear that if S contains the edge x then |Fg _1 | = 0, since 
we have a row of all zeros, and similarly for \Fg\. Since we are only interested in the 
case where the product is non-zero we can therefore just consider the the case where 
S does not contain either x or y. By the above the matrix |Fs Uy | = ±1 if T /(S U y) 
is a tree, and zero otherwise. A minor expansion in the row corresponding to y has a 
single entry of 1 (in the (c, c) position), and the rest of the entries are zero. The minor 
corresponding to the sole nonzero entry is Fg , giving \Fg\ = |Fs Uy | = ±1. Similarly 
we have that | 1 1 = — |Fs Ux | = ±1, the minus sign arising since the nonzero entry 

is in the (c— l,c) position. 

To evaluate the relative sign between |FJ| and \Fg\ we give a series of elementary 
row operations to reduce one to the other. We let S be an edge subset as above and 
G be the c x (c + 1) matrix consisting of the columns of F from S U y U x. Since 
G has a c x c submatrix of full rank it has a one dimensional kernel, which can be 
expressed naturally in terms of the cut-space: Removing the edges SUiUj/ from 
the graph gives a graph with two components, denoted A and B. For each edge in 
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SUxUy we associate the following vector w: if the edge points from A to B we assign 
the corresponding entry of w to be +1, if the edge points from B to A we assign — 1, 
otherwise we assign the entry 0. This vector must be in the kernel of G, since each 
entry of Gw gives the number of times that the cycle leaves A minus the number of 
times that it leaves B. This allows us to express the c + 1 column of G in terms of 
the first c columns of G, and thus gives a sequence of elementary row operations to 
reduce | FfT 1 1 to \Fg\. The only one which influences the sign of the determinant G 
is the orientation of the edges x and y relative to the components A and B : if the 
edges x and y both point from A to B, or vice verse, then FfT 1 1 = — \F§\, if one edge 
points from A to B and the other B to A thenlEg -1 ! = |FJ|. Equivalently T/S has 
a cycle containing both edges x and y. If this cycle traverses both edges in the same 
sense the determinants have the same sign, otherwise they have opposite signs. 

□ 

Example 3.13. One graph for which it is straightforward to compute everything 
is the complete graph, K n . We consider the case where the graph has two red edges 
(x and y) and the rest of the edges have weight 1. In this case it is clear that there 
are two topologically distinct situations, the case where the two edges share a vertex 
and the case where they do not. In the case where the two red edges share a vertex it 
is straightforward though somewhat tedious to check that the polynomial is given by 

Pk u (x, y) = 3n n_4 a;y + (2n — 3 )n n ~ 4 x + (2 n — 3 )n n ~ A y + (n — l)(n — 3 )n n ~ A 

and the discriminant is then A K n = n 2n ~ 6 . 

In the second case, where the edges do not share a vertex, the polynomial is given 
by 

P Kn (x, y) = 4n n ~ 4 xy + (2 n - 4)n n - 4 z + (2n - A)n n ~ A y + (n - 2 ) 2 n n ~ 4 

and the corresponding discriminant vanishes: A k„ =0. (Of course, this vanishing is 
implied by Corollary \3.1\ there are many automorphisms of the complete graph which 
flip one red edge and leave the other fixed.) 

4. Multiple eigenvalue crossings — R > 2. Let us now consider the case 
where there are R red edges, R > 2. As above, the crossing polynomial |2.7| is of the 
form 

A4(T) = £ Si(T)x t . 
ie 2 R 

We will find it convenient to index subsets of 2 R by the corresponding binary sequence, 
and we will write A/ = <S/(T) for brevity. 

The variety Z(T) := {x £ M(T){x) = 0} can be reducible in many dif¬ 

ferent combinations, but here we concern ourselves only with the maximal notion 
of reducibility, i.e. we consider the case when Z(T) can be written as the union of 
hyperplanes and the crossing polynomial can be written as a product of linear factors. 

It is clear that this is a very non-generic case: a general polynomial M(T) is 
determined by 2 R different coefficients, whereas an M(T) that factors is determined 
by R+ 1 different coefficients. This suggests that a necessary and sufficient condition 
for M(T) to factor should be the vanishing of 2 R — f?— 1 functions of the coefficients. 
The purpose of this section is to prove this result, and give a nice description of the 
2 r — R — 1 discriminants whose vanishing is equivalent to factorization of M(T). 

Definition 4.1. A wildcard sequence is any binary sequence on 2 R with any 
two of the digits replaced by © ’s. Let w be such a wildcard, and it corresponds with a 
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set of four binary sequences in 2 R in the canonical way: simply replace the two ©’s 
by the four sequences {00,01,10,11}, call this set S w C 2 R . For any wildcard w, we 
define 


P w (x) = ^2 A b x b . 
bes m 

We also define A w , the discriminant of P w , in the obvious manner: it is the product 
of the coefficients of the highest and lowest order terms, minus the product of the 
coefficients of the terms of middle order. 

Example 4.2. If R = 5 and w = 00©1©, then 

Pwi x ) = AoOIO^T + A00011X4X5 + ^00110^3^4 + 


and 


A w — AooioAom — Aoon Aono- 


The main result of this section is that the crossing polynomial factors into a 
product of linear factors if and only if “enough” of the reduced discriminants are zero, 
and this corresponds to the combinatorial computations of the previous section in a 
straightforward manner. We first show one direction of this theorem: 

Proposition 4.3. //A4(T) factors into R linear factors, then A w = 0 for any 
wildcard w. 

Proof. We assume that the crossing polynomial factors, so we factor it and then 
pull out all of the constants in the following manner: 

R 

A4(r)(a;) = aY\[CiXi + 1). 

i—l 

We think of C = (C\,... Cr) as a vector, and given any vector p = (pi,.. ■ ,Pr), we 
write 


S~ip S^P 1 S~lP2 S'iPr 

O — Ojl ‘ 0 2 ‘ ‘ ‘ • 

From this notation, we see that for any b £ 2 R , A b = aC b . 

Let w be a wildcard, and denote the positions of the ©’s by i < j. Let b be the 
binary sequence corresponding to w where we replace both ©’s by 0’s, and then the 
four binary sequences correspond to w are {b, b + e*, b + ej , b + e, + e^}, where e* are 
the standard basis vectors. Then 

A w = c b+Ci+e iC b - c b+ei C b+ej 


and this is clearly zero. □ 

Remark 4.4. The previous proposition shows that a complete factorization of 
M (T) implies that ( R )2 R ^ 2 discriminants vanish. While it is not hard to see that 
the converse is true, we can prove something stronger. We expect that not all of 
these discriminants are independent, and the naive count suggests that we need only 
2 r — R — 1 conditions in order to guarantee that M(T) factors. In the next part we 
define a subset of2 R — R—l whose vanishing guarantees that M(T) factors. 

Definition 4.5. For any R, we define Wr (a stacked deck) to be the set of all 
wildcards of length R satisfying the following properties: 
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• The two ©’s may be placed in any positions. 

• All of the bits before the second © must be zero. 

• The bits after the second © are unconstrained: they may be either 0 or 1. 
Proposition 4.6. The set Wr has 2 R — R — 1 elements. 

Proof. Let in be a wildcard where the second © is in slot j. Then the first © can 
be chosen in any of j — 1 positions, and all of the entries after j are free. Thus there 
are 2 R ~i(R — 1) wildcards where the second © is in the jth slot, and therefore 

R 

\W R \ = ~ ty R ~ j =2 r -R-1 . 

3 =2 


□ 

Example 4.7. For R = 3, we have \W 3 \ = 8 — 3—1 = 4. The elements of W 3 
are 

W 3 = {©©1, ©©0, ©0©, 0©©}. (4.1) 

For R = 4, we have \W 4 \ = 11. The elements ofW 4 are 

W 4 = {©©11, ©©10, ©©01, ©©00, ©0©0, ©0©1, ©00©, 0©©0,0©©1,0©0©, 00©©} 

(4.2) 


Lemma 4.8. Every sequence in Wr is of two types: 

• it is of the form wO or w\ for some w £ Wr_\, 

• it is length R, it ends with a ©, and all non-© entries are zeros. 

Proof. It is clear that any sequence described above is in Wr. Then we simply 
need to count all of the entries described above. There are 2Wr-\ of the first type 
of sequence, and R — 1 of the second type. Thus we have obtained 2Wr-\ + R — 1 
sequences of length R , and this the same recursion relation obtained by the sequence 
\Wr\r=3,...- D 

Theorem 4.9. IfT is a weighted graph with R negative edges, then a necessary 
and sufficient condition for Z(T) to decompose as a union of hyperplanes is that 
A w = 0 for all w £ Wr. 

Proof. The necessary direction of this statement has been proved already in 
Proposition |4.6[ so we consider sufficiency here. We will use a proof by induction, 
with base case R = 2; this case was established in Lemma [iPl 

Assume that the theorem is true for R—l, and that = 0 for all w £ Wr. Let 
bR C 2 r be all of those sequences with a one in the last slot, and let Pr : 2 R —> 2 R ~ 1 
be the projection on sequences that forgets the last slot. Then 

M(T) = £ S/(r)^ 7 = xr ^2 a i x P r< ' i ' > + ^2 a i x P r( ' i ' > 

Ie 2 R l£b R Igb R 


We write 


/r(z) = a i xPr{I) i 9 r( x ) = a i xFr{I \ ( 4 - 3 ) 

I^R I&bR 


giving 


A4(r) = x R f R ( x) + g R (x). 
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Note that /r(x), gn(x) are functions only of X\,... ,xr-\. 

We can write fn(x) as in (4.3), but we can also write this as 


fn( x ) = 51 Anx 1 - 

/ e 2«- 1 

Note then that the condition that /r{x) fully factor is then that A wl = 0 for all 
w £ Wr- i, but by Lemma 4.8 wl £ Wr for all w £ Wr-\. The argument for gR(x) 


is similar using wO for w £ Wr- i. Therefore, f r{x ) and gR{x) fully factor into R— 1 
linear terms. 

Now, write v,w as the vector of coefficients of fR,gR, respectively. Each of these 
vectors are of length 2 R ~ 1 . We have shown that each corresponds to a factorizing 
polynomial, and we needed | Wr-i | = 2 R ~ 1 — R conditions to establish this. Thus v, w 
each only have R degrees of freedom. Therefore, exactly R — 1 additional conditions 
guarantee that v, w are linearly dependent. Note that there are R— 1 sequences of the 
second type in Lemma |4.8[ and moreover, notice that each of them involve a pair of 
coefficients that does not appear in any of the others, e.g. the sequence with its first 
© in the ith slot involves the coefficients of XiXR,Xi. Thus they are all independent, 
and these conditions are sufficient to guarantee the collinearity of v and w. 

We have shown that Jr(x), gR(x) each fully factor into R — 1 linear terms, and 
that they are scalar multiples of each other (say gR.(x) = j3fR(x)). Then we can write 

M(T) = (x R + /3)f R (x), 

and this is clearly a product of R linear factors. □ 

Finally we note that the arguments of the previous section tell us how to interpret 
& as a signed count of spanning forests in a certain derived graph, and this is proved 


the same way as Theorem 3.6 
Proposition 4.10. LetT 
way: for each 1 digit of b the corresponding edge of T is contracted, while for each 0 
digit of b the corresponding edge of T is deleted. Then 


b be the (multi)graph derived from T in the following 




, FSF< 2 ) 


(- 1 ) 


n(F) 


5. Stability Estimates. Here we prove a sufficient condition for stability of a 
Laplacian in terms of the “worst edge”. 

Theorem 5.1. Let T be a weighted graph with fixed positive weights with T + 
connected and R red edges. For each red edge e*, let T^ = T{q be as defined as in 
Definition\2.1C\ Define 


Ui = 


Wf(r,) 

A4(T + )' 


Then, for any t = (ti, t %,..., t#) with 


||t|L < min u>i , 
i 

define Gt as the graph where we associate weight —ti to edge has spectral index. 
Then the spectral index of C(Gt) is (n_, no, n+) = (N — 1,1,0). 
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Proof. For each i, let e, be the standard basis vector with a 1 in slot i. If t = ae,, 
then for all a < t ,, L(Gt) is negative semidefinite with N — 1 negative eigenvalues. 
More generally, if ||t||^i < min,; u>i , then t is a convex linear combination of the 
basis vectors described above, and is thus a convex linear combination of positive 
semidefinite operators, and thus is itself positive semidefinite. Moreover, if we restrict 
to the subspace l - , then the previous sentence is still true after replacing “positive 
semidefinite” with “positive”. □ 

Remark 5.2. Basically, this just tells us that if we can rank the red edges from 
“best” to “worst”, then as long as we would not lose stability by putting all of the 
weight on the worst edge, then we can redistribute the same amount of negative weights 
however we choose and still retain stability. 





number of edges 


Fig. 6.1. This figure contains several observables for Gm(N, M) versus M, where we choose 
N = 10 in the top row and N = 50 in the bottom. In the left figure, we plot both the probability that 
A = 0 (black) and the probability that r_|_ is not connected (red). On the right, we plot the mean 
and standard deviation of the gap conditioned on T+ being connected and A 7= 0. 


6. Numerical Results. In this section we present a collection of numerical 
results for random matrices. We are particular interested the probability distribution 
of the nondimensionalized “gap” of a random matrix with two red edges. In all of the 
simulations performed, we have used the Erdos-Renyi random graphs G m(N, M), the 
uniform random graph with N vertices and M edges. (The specific definition for the 
distribution on these graphs is as follows. Fix N, order the N(N — l)/2 edges in some 
manner, and then choose M of them without replacement.) To make this a signed 
graph, we then uniformly choose two of these edges to be red. (We also performed 
all of the simulations presented here for the graph model Gp(N,p) (in this case, one 
chooses each edge to be present independently with probability p) and the results are 
almost indistinguishable. In the interests of space we do not present these here.) 

We see a variety of interesting behaviors in these numerics. First, we plot several 
coarse observables of this ensemble in Figure 6.1 The first quantity of interest is the 
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Fig. 6.2. Histograms of the size of the gap for the Gm(N, M) model. 



Fig. 6.3. This figure is an overlay of several conditional distributions for N = 10, M = 30, 
with specifics described in the text. In particular, the histogram in the top left frame of Figure \6.2\ 
is a convex linear combination of the various histograms presented here. 


probability that the gap is zero, and the probability that it is undefined. Recall that 
the gap will be undefined whenever the graph is disconnected, since in this case An 
is zero (q.v. (3.2)). Of course, this is more likely for a small number of edges, and 
we see these curves decreasing monotonically. Interestingly enough, the probability 
of a zero gap is actually not monotone, but in fact turns around: it is most likely 
when there are very few edges and when there are many edges. In fact, we see that 
for sparse graphs, the vast majority of the connected graphs have zero gaps. In all 
of these pictures, the number of edges increases until we reach the complete graph; 
as proved in the text, the event of having a zero gap for the complete graph is the 
same as having the two red edges not share a vertex (q.v. Example 3.13). A simple 
combinatorial argument shows that this probability is 1 — 4/IV, giving 0.6, 0.92 in the 
two graphs plotted in Figure [6(1] 

In Figure |6.2[ we plot the histograms of the gap in a variety of ways. As in 
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Figure |6.1[ the top row corresponds to N = 10 and the bottom row to N = 50. In 
the left column, we have plotted the histogram for a single value of M; in the top 
left, we have N = 10, M = 30 and in the bottom left we have N = 50, M = 605. 
In the right column, we plot a heatmap of the histograms for various M; in the top 
right, we plot all M in the range [10,45] and in the bottom right, we have M going 
from 55 to 1225 by jumps of 10. Note that in each row, the left frame is a horizontal 
“slice” of the right frame. Each “slice” in the right figure contains 10 5 realizations 
of random graphs; for the specific histograms in the left column we simulated 10 7 
random graphs. 

What is most striking about this figure is the lack of normality of the distribution, 
yet, at the same time, the clear impression that they are multimode normal. Moreover, 
the heatmaps certainly suggest that these modes are coherent as a function of M. 
Again, we know that the distribution for the complete graph is degenerate in that it 
can take on only two values: zero (if the edges are disjoint) and N 2N ~ 6 (if the edges 
share a vertex). Since we are plotting the logarithm of the gap, the zero is represented 
by a mass at — oo (which is here binned at —10). So we see that only one nonzero 
mode “survives” as M goes to the complete graph, and this is the one corresponding 
to adjacent edges. 

It is natural to ask what these modes correspond to, and we have decomposed the 
distribution in Figure 645 We have done this for the case N = 10, M = 30, so that 
this is a decomposition of the histogram in the top left frame of Figure [672] What we 


have done is as follows: for all of the random graphs generated for N = 10, M = 30, 
we condition these in various ways. Those where the red edges are adjacent we put in 
class “adj”. For those in which the red graphs were not, we computed four pairwise 
path distances between the vertices of the two red edges. To be more specific, if the 
two red edges were x = (x±,X 2 ) and y = ( yi,V 2 ), we computed dG + (xi,yj), i,j = 1,2 
that is to say, we computed the graph distance between the nodes using only black 
edges. If, for example, all of these vertices were adjacent in the black graph, we put 
it in class “1111”; if, for example, three of them were adjacent, but one pair had a 
G .|_ distance of 2, we put it in class “1112”, etc. We then plotted the conditional 
distributions for each class, and we see the conditional distributions are close to log¬ 
normal (recall that the horizontal axis is always the logarithm of the gap). 

In summary, while the distribution of all graphs with a fixed number of edges is 
multimodal, if we consider those graphs with a fixed number of edges M. and two 
red edges, and we condition on the geometry of these red edges, then the conditional 
distributions are each quite close to log-normal. 
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